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ABSTRACT 

We argue that the lack of observed damped Lyman a (DLA) systems that simultaneously have 
high H I columns densities and high metallicities results naturally from the formation of molecules in 
the cold phase of a two-phase atomic medium in pressure balance. Our result applies equally well in 
diffuse systems where the ultraviolet radiation field is dominated by the extragalactic background and 
in dense star-forming ones, where the local radiation field is likely to be orders of magnitude higher. 
We point out that such a radiation-insensitive model is required to explain the absence of high column 
- high mctallicity systems among DLAs observed using gamma-ray burst afterglows, since these are 
likely subjected to strong radiation fields created by active star formation in the GRB host galaxy. 
Moreover, we show that the observed relationship between the maximum atomic gas column in DLAs 
sets a firm upper limit on the fraction of the mass in these systems that can be in the warm, diffuse 
phase. Finally, we argue that our result explains the observed lack of in situ star formation in DLA 
systems. 

Subject headings: galaxies: ISM — gamma rays: bursts — ISM: molecules — ISM: structure - 
quasars: absorption lines — stars: formation 



1. INTRODUCTION 

Damped Lyman a (DLA) systems arc clouds of neutral 
atomic hydrogen with column densities iV(H i) > 2 x 10 20 
cm~ 2 that are detected as absorbers against bright back- 
round quasars (QSOs) or gamma-ray bursts ( GRBs ) 
Wolfe et all [20051 : iProchaska et~aTl I2005L 120071 . 120081 ). 
These systems comprise the bulk of the neutral gas in 
the universe at redshifts up to at least z ~ 5. Because of 
their ubiquity, the study of DLAs provides vital clues to 
the distribution of gas and metals, and potentially also 
star formation, in the universe. 

Observations of DLAs show a clear zone of exclusion: 
none are observed with both high metallicity and high 
H I column density. One possible explanation for this ef- 
fect is that lines of sight with large column densities and 
metallicities produce large dust extinctions that might 
lead to exclusion of the background QSO from optically- 
select ed samples dBoisse et al.lll998l iPrantzos fe Boissien 
2000). However, statistical analysis of the optically- 
selected QSO-DLA sampl e suggests that few DLA s are 
missed due to extinction (jPontzen fc Pettinjl l2009f ) . and 
radio-selected QSO-DLA samples d o not differ signifi- 
cantly from optically-selected ones (lEllison et all I2001L 
120051 lAkerman et"aLll2005l : Uorgenson et al.ll2006l ). 

An alternative hypothesis to explain the zone of ex- 
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elusion is that above some threshold total hydrogen col- 
umn density, which decreases with increasing metallicity, 
gas forms molec ular hydrogen that is not detectable in 
Ly a absorption jS chavc 200 It iHirashita fe Ferrarall2005l : 
iHirashita et al.l l2006f ). Large amounts of molecular gas 
are not observed in these DLAs because molecular clouds 
have a very small covering fraction and are unlikely to 
be se en along random sightlines (|Zwaan fc Prochaskal 
120061) . although in some cases trace amounts of molec- 
ular gas have been de tected (e.g. iLedoux et al.1 120031 : 
iNoterdaeme et al.l 120081 ) . In one GRB-DLA a substan- 
tial column of molecu lar hydrogen has been detected 
(|Prochaska et a l. 2009), and, as we show below, this sys- 
tem is unique among DLAs in its high column density 
and mctallicity. 

While the molecule-formation hypothesis avoids the 
problems of the dust bias explanation, it also has sig- 
nificant weaknesses. In DLAs we can observe only col- 
umn density and metallicity, so previous authors have 
been forced to assume values, which may be incorrect, 
for other quantities such as total gas volume density 
and radiation field that influence the molecule frac- 
tion. For example, ISchavd (|2001f ) assumes a Lyman- 
Werner (LW) rad i ation field within a factor of 3 of the 
lHaardt fc Madaul (|2001f ) UV background at z = 3, which 
is 50 times smaller than the Solar neig hborhood value, 
while the weakest radiation field that fHirashit a et al.l 
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(|2006f ) consider is 2 times larger than in the Solar neigh- 
borhood. Observationally-infcrred ra diation fields in 
DLAs span this full ra nge and more ([Tumlinson et all 
l2007tlWolfe et al.ll2008f ). so neither assumed value works 
for the full DLA population. 

Our goal in t his paper is to explain t he zon e of exclu- 
sion using the IKrumholz etall (|2008l l2009al . hereafter 
KMT08 and KMT09) theory of the atomic to molecular 
transition in galaxies, which does not depend on unob- 
servable quantities such as the gas v olume density and 
UV radiation field. Instead, IKMT09I show that to good 
approximation the molecular fraction in a cloud depends 
only on its column density and a single dimensionless 
parameter, which combines the volume density, radiation 
intensity, H2 formation rate coefficient, and dust opacity, 
and that the two-phase nature of the atomic interstellar 
medium imposes strong constraints on the values this 
parameter can take. This model explains the observed 
molecula r fractions and star formation r ates in nearby 
galaxies (|KMT09l : IKrumholz et aTll2009bf ). and in §[2] we 
apply it to DLAs. In §[3] we demonstrate that the ob- 
served zone of exclusion sets strong constraints on the 
fraction of warm atomic gas in DLAs. Finally, in §[4] we 
discuss the implications of our work. 

2. MOLECULE FORMATION AND THE ZONE OF 
EXCLUSION 

2.1. Theoretical Model 0/H2 Formation 

We refer readers to IKMT08I and IKMT09I for a full 
derivation of the formalism, and here simply summarize. 
Consider a spherical cloud of gas immersed in a uniform, 
isotropic dissociating radiation field. The outer parts 
are kept predominantly atomic by the radiation, but as 
one moves inward dissociating photons are absorbed by 
H2 molecules and dust grains. At some depth into the 
cloud all dissociating photons have been absorbed, and 
there is a sharp transition to predominantly molecular 
material. The fraction of the cloud radius at which this 
transition occurs depends on two dimensionless numbers: 
Td = Na ( i, the dust optical depth of the cloud, and 
X = /dissCdC-Eg / ' (nmTV) , the dimensionless strength of 
the radiation field. Here N is the center-to-edge col- 
umn density of the cloud including atomic and molecu- 
lar material, <7d is the dust cross section per H nucleus 
to photons in the dissociating LW bands, /diss ~ 0.1 is 
a quantum-mechanical constant describing the approx- 
imate probability of dissociation per absorption, Eq is 
the number density of photons in the LW bands of the 
background dissociating radiation field, nni is the num- 
ber density of gas in the atomic envelope of the cloud, 
and 1Z is the rate coefficient (in units of length 3 time -1 ) 
describing formation of hydrogen molecules on the sur- 
faces of dust grains. In the limit N — » 00, the Hi col- 
umn density JV(Hi) approaches a finite maximum value 
that depends only on x an d on the dust cross section 
Od (which in turn depends on metallicity) in the atomic 
shielding region around the cloud. In effect, a fixed op- 
tical depth of material can absorb the entire LW photon 
flux, so any additional gas is molecular. This produces 
the observed zone of exclusion: a metallicity-dependent 
maximum Hi column. 

IKMT09I point out that x cannot vary strongly be- 
tween galaxies. Both o~d and 1Z measure the total sur- 
face area of dust grains in the gas, so u^/'R, is insensi- 



tive to changes in dust abundance or size distribution, 
and thus varies little with environment. Moreover, in 
a two-phase atomic medium, the H 1 density nm that 
determines when molecules form is the density of the 
cold neutral phase, tt-cnm- This is because the effec- 
tive LW opacity of a fluid element provided by H2 ab- 
sorption is proportional to its density (cf. equation 8 of 
IKMT08T ). Thus the low density of the warm phase guar- 
antees that it provides negligible self-shielding compared 
to the cold gas, and we care almost exclusively about 
ncNM- Of course the WNM does contain dust, but gas 
only becomes predominantly molecular once the ambi- 
ent UV radiation field has been attenuated by a fac- 
tor of ~ 10 3 — 10 4 . This level of dust attenuation re- 
quires a color excess E(B — V ) > 0.6, larger than the 
highest known E (B — V) (e.g. Uunkkarinen et all 120041 : 
I Wild et aLll2006l) and ~ 2 orders magnitude above the 
mean (lEllison et al.|[2005l : IVladilo et"at| [2008). Thus we 
can safely neglect the contribution of WNM dust shield- 
ing in favor of CNM self-shielding. The ratio Eq/ticnm 
is tightly constrained by the thermodynamics of the gas 
and the requirement of pres sure balance between the two 
phases (jWolfire et a l. 2003): a reasonable approximation 
is E'Jncau ~ (1 + 3.1Z'°- 365 )/93 cm 3 , where Z' and 
E' are the metallicity and FUV radiation intensity Eq 
normalized to their values in the Solar neighborhood. 1 
Together, with the invariance of o~d/Tt, this (weak) de- 
pendence of Eq I hcnm on Z' gives a dimensionless radia- 
tion intensity x ~ 0.77(1 + 3.1Z' 365 ) that depends only 
on the metallicity of the gas. 

The implication of this result is that, in the CNM, the 
dimensionless radiation strength x does not depend on 
the absolute FUV radiation field. Any change in radia- 
tion intensity induces a countervailing change in density. 
This is why there is a zone of exclusion for all DLAs 
despite the huge range in radiation intensities inferred 
within them. To calculate this effect quantitatively, if 
we assume that dust opacity o& scales with metallicity, 
then in the KMT formalism the molecular mass fraction 
is approximately given by 
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where N c is the column density of cold gas (i.e. in- 
cluding CNM and molecular gas, but excluding warm 
atomic gas), s = ln(l + 0.6x)/(0.045A^20^')i ^20 = 
N c /[10 20 H nuclei cm" 2 ] and S = 0.0712(0.1s- 1 + 
0.675)~ 28 . ( This appr oximation is slightly different than 
that given in IKMT09I : the two agree closely for clouds 
that are substantially molecular, but this one is more ac- 
curate at low /h 2 - McKee & Krumholz, 2009, in prepa- 
ration.) This expression gives the mass of the spherical 
molecular core of the cloud, which is surrounded by a 
shell of atomic gas whose density is lowe r than th at of 
the molecular gas by a factor <p mo \ , which IKMT09I show 
is typically ~ 10. The covering fraction of the molecular 
sphere is 
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1 Following Drainc (1978), wc take the LW radiation intensity in 
the Solar neighborhood to have a value that produces a fr ee-space 
dissoc iation rate of 5.43 X 10 -11 s _1 ; this is 1.6 times the lHabind 
il968h field. 
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Fig. 1. — Hi column density TV (Hi) versus normalized 
metallicity Z' and zinc abundance log(Zn/H) + 12, computed 
for cold gas. We show lines of constant H2 covering fraction 
(solid blue, ch 2 increasing with ./V(Hl)), lines of constant color 
excess (dott ed green, E(B — V ) incre asing with iV(Hl)), QSO- 
DLAs from IHerbert-Fort et al.l j2006ft , Kaplan et al. (2009, in 
preparation), and Dessauges-Zavadsky et al. (2009, in preparation) 
(black circles), GRB-D LAs without H2 detections (red diamonds; 
IProchaska et al.ll2007ft . and the G RB080607-DLA with a n H 2 de- 
tection (purple diamonds with line, IPro chaska et al. 20091). For the 
QSO-DLAs, filled circles indicate detections of metals and open cir- 
cles indicate the la upper limits on metallicity. See the main text 
for discussion. 

Obviously a spherical ball is a great oversimplification of 
the complex geometries of atomic-molecular complexes, 
but ch 2 is a useful general indicator of the fraction of the 
area that is likely to be covered by molecular material. 

For a given metallicity Z' it is trivial to numerically 
invert equation to calculate the total cold gas col- 
umn density N c for which the molecular covering fraction 
reaches a particular value ch 2 - The mean atomic col- 
umn density is then iV(Hl) = (1 — fn 2 )N c . This defines 
a locus of points in the iV(Hl), Z'-plane corresponding 
to the specified ch 2 ■ The maximum H I column density 
corresponds to the limit ch 2 — > 1, because this corre- 
sponds to an infinite slab illuminated by the external 
radiation field. At this point we must mention two im- 
portant caveats. One is that we assume that the atomic 
gas in DLAs is in two-phase equilibrium, which may not 
be true for all of them. The other is that this method 
allows us to constrain only the cold H I column density. 
In principle much larger warm gas column densities are 
possible (§ [3j) . 

2.2. Comparison of Models and Observations 

In Figure[l] we plot our derived values iV(H i) versus Z' 
for molecular covering fractions from ch 2 = 0.01 — 1. 2 Wc 
also show line s of constant E{B — V), computed using a 
iDraind (|2003T ) Ry = 3.1 extinction curve scaled by metal- 
licity, giving E(B-V)/N(Ri) = 1.65 x 10" 22 Z' cm 2 and 
A v /N(Ri) = 5.32 x 10~ 22 Z' mag c m 2 . We compare to 
observ e d QSO- and GRB - DLAs f romlHerbert-Fort etal] 
(|2006| ). IProchaska et~aTl (|2007l 120091), Kaplan et al. 
(2009, in preparation), and Dessauges-Zavadsky et al. 
(2009, in preparation). For the Dessauges-Zavadsky 
et al. sample we derive metallicities from zinc abun- 
dance: logZ' = [Zn/H] = log(Zn/H) - log(Zn/H) Q , 

2 For ch 2 = 1 we use the KMT formalism to solve for /h 2 
numerically rather than using equation JTJ. This is necessary be- 
cause as en — > 1, 7V(Hl) — > (1 — fn 2 )N c depends on depends on 
dfn 2 /dN c . Equation JT} does not give precisely the correct limit 
for this quantity. 
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Fig. 2. — Same as Figure [T] plus blue lines showing the max- 
imum values of N(Hl) and Z' assuming that warm gas fractions 
/ w = 0.0 (solid), / w = 0.5 (dashed), and / w = 0.9 (dot-dashed). 
QSO-DLA points shown in gree n are those with measured spin tem- 
peratures {Kanckar et al. 20091). The open purple dia mond shows 
the estimated GRB080607-DLA cold Hi column (see § |2~3)) . 

where log(Zn/H) + 12 = 4.63 (|Loddersl |200~1 . For all 
other data we use the metallicity reported by the authors. 
To avoid possible issues arising from cither ionization cor- 
rection or metallicity evolution with redshift, we exclude 
DLAs with logN(Hi) < 20 and redshift z < 1.7. 

As the Figure shows, the zone at high tV(Hi) and 
Z' where no DLAs lie (except that associated with 
GRB080607, which wc discuss below), corresponds well 
to the predicted zone of exclusion. The molecular cover- 
ing fraction declines sharply away from the ch 2 = 1 line, 
so all DLAs but GRB080607 lie below c H ., = 0.06 This i s 
consistent with the results of lZwaan fc Pro chaska ( 2006h . 
who conclude that detection of true molecular clouds 
in DLAs is unlikely because the molecular material has 
a small covering fraction. Trace amounts of molecu- 
lar hydrogen have been discovered in some DLAs (e.g. 
iLedoux et all 120031 iNoterdaeme et"aT1l2008| ). but these 
low molecular columns almost certainly correspond to 
H2 spatially mixed with cold atomic gas, rather than true 
molecular clouds. The KMT formalism approximates the 
atomic-molecular transition as sharp, so it does not ap- 
ply to these systems. We defer discussion of them to 
future work. Also note that the observed distribution 
falls off sharply at logZ' > 0, and at logA^(Hl) > 22 
independent of Z'. Molecule formation cannot explain 
these features. 

2.3. The DLA Associated with GRB080607 

The DLA asso ciated with GRB080607 
()Prochaska et all I2009D . the only DLA inside the 
zone of exclusion, is also the only DLA to show signif- 
icant columns of H2 and CO. We plot this detection 
at two metallicities derived in different ways. The 
\ogZ' = —0.2 point corresponds to the oxygen abun- 
dance [O/H]. The logZ' = —0.94 point is derived using 
the observed visual extinction Ay ~ 3.2 and Hi and H2 
columns logiV(Hi) = 22.7 and \ogN(R 2 ) = 21.2, and 
adopting the same metallicity-dependent Ay /N(Hi) 
ratio as in § 12.21 Since the molecular fraction depends 
on solids that can catalyze H2 formation and absorb LW 
photons, the latter estimate is probably the relevant one. 
The solid content of the atomic gas may be even lower 
if a disproportionate share of the observed extinction 
comes from the molecular material. 

While the detection of molecules is consistent with the 
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DLA's presence in the zone of exclusion, we have not yet 
explained why its column is only 6% molecular. Dissoci- 
ation by the GRB afterglow is unlikely to be the explana- 
tion. The hard afterglow spectrum would produce nearly 
coincident ionization and dissocia tion fronts with littl e 
atomic hydrogen between them (jDraine fc Had l20()3) . 
and excited H 2 just outside the dissociation front would 
produce stron g absorption features i n the DLA that are 
not observed ([Prochaska et al.|[2009t however, see Shef- 
fer et al., 2009, in preparation). Observations instead 
suggest that the molecular cloud is at least 100 pc from 
the GRB, so our line of sight must pass through the host 
galaxy's disk at a glancing angle. Thus a majority of the 
material along the line of sight should be warm atomic 
gas that is unrelated to the molecular cloud, and, as dis- 
cussed in § 12. U provides no shielding to it. 

Let the DLA sightline contain warm and cold (atomic 
plus molecular) material with column densities _/V w and 
N c . The molecular column will be 

l.W(H 2 ) = f H2 (N c ,Z')N c , (3) 
where the factor of 1.8 accounts for the difference in mean 
number of particles per unit mass between atomic and 
molecular gas. For logZ' = -0.94 and log A(H 2 ) = 21.2, 
solving this equation yields logA c = 22.0, which implies 
/h 2 = 1.8Nu 2 /N c = 0.31 and logA w = 22.6. Thus we 
can explain the observed Hi and H 2 column densities if 
18% of the H I column consists of cold gas and the remain- 
ing 82% is warm. The open purple diamond in Figure 
[5] shows where this system would fall in the N(Hl),Z' 
plane if we counted only cold Hi. The path length 
through the warm gas is L = 13(?t-wnm/1 cm" 3 )" 1 kpc, 
where tt-wnm is the volume density. For nwNM ~ 10 -05 
cm" 3 , a typical value in the Solar neighborhood, this 
would require L ~ 40 kpc, but in a two-phase medium 
the equilibrium WNM density is close to linearly pro- 
portional t o the intensity of t he FUV radiation field in 
the galaxy (jWolfire et alj2003j) . Since the FUV radiation 
field in other GRB host galaxies is ~ 10 — 100 times larger 
than the Solar neighborhood value (jTumlinson et all 
12007ft . we expect tiwnm ~ 10 cm -3 , giving L ~ 1 kpc. 

3. CONSTRAINING THE WARM GAS FRACTION 

Warm gas adds H I column without creating molecules, 
possibly moving systems to the right into the exclusion 
zone. To illustrate this, in Figure [2] we show the same 
observed systems as in Figure [TJ compared with the line 
ch 2 = 1 calculated assuming that a fraction f w of the 
observed H I column is in the form of warm gas that pro- 
vides no shielding. These lines indicate the maximum 
values of iV(Hi) that can be observed for a given mctal- 
licity Z' for the indicated value of / w . 

DLAs with large total gas columns, high Z', and large 
/ w , would not form molecules, so they would be observed 
to the right of the f w = line. The fact that this region 
is unpopulated by QSO-DLAs suggests that the combi- 
nation of large gas column, large Z' and large / w must 
be very rare for them. It may be more common for GRB- 
DLAs, since one of six systems falls within the exclusion 
zone. Of course even if such large column-large Z' DLAs 
did exist, the distribution would likely be cut off at some 
point by dust extinction effects, and even in the absence 
of CNM a sufficiently large WNM column could provide 
enough shielding for molecules to form. A precise es- 
timate of the warm gas fraction permitted by the data 



would require a careful analysis of these effects. Nonethe- 
less, the fact that dust extinction alone cannot explain 
the zone of exclusion strongly argues that / w cannot be 
anywhere near unity for the DLAs with the highest val- 
ues of i V(Hi) and Z ' . Our conclusion is consistent with 
that of I Wolfe et "all (|2008l ). who argue based on C II* 
absorption measurements that ~ 1/2 of DLAs contain 
significant cold gas. 

However, we emphasize that our model offers no con- 
straints on the warm gas fraction in DLAs with total 
column densities and metallicities that place them below 
the / w = 1 line. Thus our resu l ts are consistent with 
the conclusion of iKanekar et al.l <|2009f ) . who find large 
values of f w for the DLAs shown by the green points in 
Figurc[2] Significantly, all of these points are at least 1.6 
dex to the left of the f w — line. 

4. SUMMARY AND DISCUSSION 

We show that the absence of high column density, high 
metallicity DLAs results from the conversion of atomic 
into m olecular gas, confirming the proposal of iSchavd 
(|2001l ). Our model is insensitive to the DLA radiation 
environment due to the way that two-phase atomic gas 
responds to variations in radiation field. The maximum 
metallicity and column density do depend on the frac- 
tion of the atomic gas that is cold and dense, and the 
observed distribution of QSO-DLAs is inconsistent with 
the existence of a significant population of high column 
density, high metallicity QSO-DLAs dominated by warm 
gas.^ 

Given our conclusion that high column density DLAs 
must host significant amounts of cold gas, o ne might ask 
why, as pointed out bv lWolfe fc CherJ ( 20061 ). most DLAs 
cannot host significant in situ star formation. The an- 
swer is that the presence of cold gas is a necessary but not 
sufficient condition for star formation. Although DLAs 
populate the N(Hi),Z' plane up to the point where 
the molecule fraction becomes large, the vast majority 
of them are found at much lower column densities and 
metallicities, where they are not expected to have any 
significant amount of molecular gas. If stars form exclu- 
sively in molecular gas, as numerous observations now 
seem to suggest, then the vast majority of DLA columns 
should be inert as far as star formation is concerned. 

This is not to say that DLA systems do not host any 
star formation. Indeed, indirect measures of the radia- 
tion fields in some QSO-DLAs show strong evidence for 
the presen ce of a local heat so urce that is likely to be star 
formation (|Wolfe et al.ll2008f) . and in GRB-DLA systems 
there is obviously evidence for ongoing star formation. 
Our result simply suggests that the star formation must 
be taking place in other parts of these galaxies, which 
have significantly higher column densities and molecular 
contents than the sightlines we most commonly observe 
as DLAs. 
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